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ABSTRACT 

We studied the brightness and spectral evolution of the young eruptive star 
V1647 Ori during its recent outburst in the period 2004 February - 2006 Sep. We 
performed a photometric follow-up in the bands V , Rc, Ic^ H, Kg as well as visible 
and near-IR spectroscopy. The main results derived from combining our data with those 
published by other authors are as follows: The brightness of V1647 Ori stayed more than 
4 mag above the pre-outburst level until 2005 October when it started a rapid fading. In 
the high state we found a periodic component in the optical light curves with a period 
of 56 days. The delay between variations of the star and variations in the brightness of 
clump of nearby nebulosity corresponds to an angle of 61° it 14° between the axis of the 
nebula and the line of sight. The overall appearance of the infrared and optical spectra 
did not change in the period March 2004 - March 2005, though a steady decrease of 
HI emission line fluxes could be observed. In 2006 May, in the quiescent phase, the 
Hel 1.083 iim line was observed in emission, contrary to its deep blueshifted absorption 
observed during the outburst. The J — H and H — Kg color maps of the infrared nebula 
reveal an envelope around the star whose largest extension is about 18"(0.03pc). The 
color distribution of the infrared nebula suggests reddening of the scattered light inside 
a thick circumstellar disk. Comparison of the Kg and Ha images of McNeil's Nebula, 
the conical nebulosity illuminated by V1647 Ori, shows that HH 22A, the Spitzer in- 
frared source and the bright clump C of the nebula may be unrelated objects. We show 
that the observed properties of V1647 Ori could be interpreted in the framework of 
the thermal instability models of Bell et al. (1995). V1647 Ori might belong to a new 
class of young eruptive stars, defined by relatively short timescales, recurrent outbursts, 
modest increase in bolometric luminosity and accretion rate, and an evolutionary state 
earlier than that of typical EXors. 

Subject headings: stars: formation-stars: circumstellar matter-stars: individual: V1647 
Ori - stars: pre-main sequence-ISM: individual objects: McNeil's Nebula 
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Introduction 



Eruptive young stellar objects form a small, but spectacular class of pre-main sequence stars. 
Traditionally they are divided into two groups. FU Orionis-type stars (FUors) are characterized by 
an initial brightening of ~5 mag during several months or years, followed by a fading phase of up to 
several decades or a century. Their spectral type is F-G giant according to the optical spectrum, 
and K-M giant/supergiant at near- infrared wavelengths. The second group, called EX Lupi-type 
stars (EXors), belongs to the TTauri class. Their recurrent outbursts are relatively short, lasting 
from some weeks to months, and the time between the eruptions ranges from months to years. 
Their spectral type is K or M dwarf. 

The recent outburst of the young star V1647 Ori offered a rare opportunity to study phenomena 
accompanying the eruption of a low-rn ass pre-main sequence star. Following a 3-month brightening 
of A/c ^ 4.5 mag (jBricerio et al.ll2004l ). V1647 Ori reached peak light iii January 2004, illuminating 
also a conical reflection nebulosity called McNeil's Nebula (|McNeilll2004l). The obiect st ayed in high 



state for about 2 years, then started a rapid fadin g in October 2005 (IKospal et al 
exhibited already a similar eruption in 1966-67 (jAspin et al.l 120061 ). 



20051). V1647 Ori 



During the 2004-2005 outburst V1647 Ori was intensively studied from many aspects. The 



bright n ess evolution at optical and n ear-infrared (NIR) w a velengths was monitore d bv lBricefio et al 



2005 



( 2004 ) iReipurth fc Asuin I (120041 ). IWalter et al.l (120041 ). iMcGehee et al.l (120041). and lOi 



20061). The o ptical/infrared s pectr u m, and its changes were studied by 



Vacca et 



al 



la et al 
2004|) 



Walter et al] (|2004l ). lA^drews et al.l (2004l ). lRettig et al. I (l20od ). IOiha etahl (|2006l ). and loibb et al 
(I2OO6I). The structure of t he circumstellar niatter and t he physics of the out burst were investigated 



bv lMcGehee et al.l (j2004l ). Muzerolle et al 



Grosso et al 



(I2OO4I ) . andlEettig et al. I ( 20051 ) . The X-ray pr o pertie s 



2OO5I ). iKastner et aD (|2004l ). and Ixastner et al] (|2006l ) 



of the star were studied by 
These works produced a wealth of very important information on the eruption, but a general pic- 
ture is still to be worked out. Nevertheless, all authors seem to agree that: V1647 Ori is a deeply 
embedded low-mass pre-main sequence object; it is surrounded by a disk; the outburst was caused 
by increased accretion; and that the eruption was accompanied by strong stellar wind. It is still an 
open question whether V1647 Ori belongs to the FUor or the EXor class, and whether the standard 
outburst models, developed for FU Orionis events, could be applied for this case (or how should 
they be modified to make them applicable). 

In this paper we present and analyse the results of our monitoring programme on V1647 Ori. 
The data, which include optical ( V , Rc, Ic) and near-infrared {J, H , Kg) imaging, and intermediate- 
resolution optical and near-infrared spectroscopic observations, were obtained between February 
2004 and May 2006, covering the whole outburst period. Following a description of the observa- 
tions and data reduction (Sect. [2]), in Sect. 3 we present results on the photometric evolution of 
V1647 Ori. Sects. 4 and 5 are devoted to the spectroscopic variability and the morphology of 
the nebula, respectively. In Sect. 6 we discuss the physics of the outburst, and comment on the 
FUor-or-EXor question. Our results are briefly summarized in Sect. [71 
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2. Observations and data reduction 

2.1. Optical imaging and photometry 

Optical observations of V1647 Ori were obtained on 45 nights in the period 2004 February 12 - 
2006 February 1. The measurements were performed in two observatories with three telescopes: the 
IAC-80 telescope of the Teide Observatory (Spain), and the Im Ritchey-Chretien-Coude (RCC) 
and 60/90/180 cm Schmidt telescopes of the Konkoly Observatory (Hungary). 

The IAC-80 telescope was equipped with a 1024 x 1024 Thomson CCD with a scale of O'.'4325/pixel, 
and a field of view of Y'AxJ'A. For the broad-band observations we used the IAC#72, IAC#71, 
and IAC#70 filtertQ. Typically three frames were taken with exposure times of 120-300 s/frame 
each night. Additionally, on the night of 2004 January 11 deep Ha and S[II] images were obtained 
using the narrow band filters IAC7^17 and IAC7^68. Four exposures of 1200 s were taken through 
each filter. 

At Konkoly Observatory, the 1 m RCC telescope was equipped with a Princeton Instruments 
VersArray 1300B camera. The back- illuminated 1340x1300 CCD has an image scale of O'.'306/pixel 
and a field of view of 6f8x6f6. The Schmidt telescope was equipped with a 1536x1024 Kodak 
Photometries AT 200 CCD camera. The image scale was l'.'03/pixel, providing a field of view of 
24' x 17'. On b oth telescopes, the ut ilized V(RI)c filters matched closely the standard Johnson- 



Cousins system (jHamilton et al.ll2005l ). For each filter, 3-10 frames were taken on each night, and 



the integration times varied between 120 and 600 s/frame. 

Data reduction and photometric calibration. Raw frames were dark-subtracted and flat- 
fielded using the 'imred' package within the IRAfI^ environment. Dome flat-field images were 
taken each night, and sometimes sky flats were also available. In some cases, the consecutive 
frames taken through the same filter were co-added in order to increase the signal to noise. As an 
example, Fig. [1] shows R-band optical images of V1647 Ori taken with the Im RCC telescope of 
the Konkoly Observatory at three different phases of the outburst. The narrow-band Ha and S[II] 
images were bias-subtracted, flat-fielded and combined per filter. Then a scaled R-band image, 
assumed to be a proper continuum, was subtracted from them in order to show the pure emission 
line map. 

The instrumental V(RI)c m agnitudes of V 1647 Ori were determined by PSF-photometry using 



the 'daophot' package in IRAF (IStetson 1119871 ). The field of view of our telescopes allowed us to 



include several field stars to define the point spread function of the images. The use of PSF- 



^The transmission curves are available at http://www.iac.es/telescopes/tcs/filtros-eng.htm 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for the Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



http://iraf.noao.edu/ 
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photometry was important to exclude the contribution of the nebula, which becomes dominant 
especially in the V-band. In order to assure the consistency of the results obtained with the three 
different telescopes we treated our data in a homogeneous manner: the position and extension of the 
sky region near V1647 Ori were the same in all images, and the preliminary aperture photometry, 
used for scaling the PSF magnitudes was obtained using equally narrow apertures (1.5 arcsec) in 
each image. 



For photometric calibration we used the secondary standards from the list of ISemkov I (|2006l ). 
Differential photometry was performed by computing the instrumental magnitude differences be- 
tween V1647 Ori and each comparison star. The differential magnitudes were transformed into the 
standard photometric system as follows. For the comparison stars, we searched for a relationship 
between Am, the differences between instrumental and absolute r nagnitudes in the V, Rg or Iq 



band and the color indices of the stars, similarly as described in ISerra-Ricart et al.l (jl999l ). We 
note that the color indices of the comparison stars covered a large enough range to bracket the 
color indices of V1647 Ori. The data points were fitted with a linear relationship, except a few 
low-quality nights, when due to large scatter, all Am values were averaged. For a given telescope, 
the derived color terms (the slope of the fit) were consistent from night to night within the formal 
uncertainties provided by the fitting procedure. Color terms computed from the 1 m RCC observa- 
tions also agreed with values determined with high precision from observations of NGC 7790 with 
the same instrument. For calibration of the V1647 Ori data we applied the color terms determined 
for the corresponding night. Results of the optical photometry of V1647 Ori are presented in Tab. [1] 
and plotted in Fig. [2j The uncertainties given in Tab. [T] were computed as the quadratic sum of 
the formal errors of the instrumental magnitudes provided by IRAF and the uncertainties of the 
standard transformation. In those cases when the frames were evaluated separately (rather than 
co-added), the standard deviation of the magnitudes from each frame was also computed, and the 
maximum of the two types of errors was adopted. 




30 20 10 -10 -20 -30 30 20 10 -10 -20 -30 30 20 10 -10 -20 -30 

OffRA (arcsec) Orf RA (arcsec) OFfRA (arcsec) 



Fig. 1. — Examples of R-band optical images of V1647 Ori taken with the Im RCC telescope of 
the Konkoly Observatory at three different phases of the outburst. North is at the top and east to 
the left. V1647 0ri is located at the southern apex of the nebula. The color table has been scaled 
to represent the same brightness in all three images. 
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Fig. 2. — Light curves of V1647 Ori in the V (dots), Rq (asterisks), Iq (upward triangles), 
J (squares), H (downward triangles), and Ks (diamonds) bands between 2004 February 16 and 
2006 May 5 , based on the data presented in Tables [1] and [2j The uncertainties are comparable 
with symbol sizes. 
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Table 1. Optical photometry of V1647 Ori 



Date 


JD^ 




V 




Rc 






Ic 


Telescope 


20040212 


048.3 














14 


33 


(0.04) 


Schmidt 


200402 20 


056.3 






17 


01 


(0 


08) 


14 


62 


(0.04) 


RCC 


200403 06 


071.3 














14 


56 


(0.04) 


RCC 


2004 04 01 


097.2 














14 


83 


(0.08) 


Schmidt 


2004 04 02 


098.3 






16 


28 


(0 


08) 


14 


51 


(0.06) 


Schmidt 


200404 24 


120.3 


18 


92 (0.06) 


16 


91 


(0 


04) 


14 


73 


(0.07) 


IAC-80 


200408 31 


248.6 


19 


00 (0.12) 


16 


86 


(0 


10) 


15 


13 


(0.15) 


RCC 


200409 01 


249.6 














14 


89 


(0.08) 


IAC-80 


200409 02 


250.6 


19 


58 (0.08) 


17 


00 


(0 


08) 








IAC-80 


2004 09 03 


251.6 














14 


74 


(0.04) 


RCC 


2004 09 06 


254.6 






17 


08 


(0 


03) 


14 


90 


(0.04) 


RCC 


200409 07 


255.6 


19 


32 (0.06) 


17 


62 


(0 


08) 


14 


95 


(0.03) 


RCC 


20041006 


284.7 


18 


91 (0.10) 


17 


24 


(0 


06) 








IAC-80 


2004 10 25 


303.6 


19 


23 (0.10) 


17 


50 


(0 


05) 


15 


21 


(0.06) 


IAC-80 


20041101 


310.8 














14 


99 


(0.04) 


IAC-80 


2004 11 21 


330.7 


18 


86 (0.08) 


17 


08 


(0 


08) 








IAC-80 


2004 11 24 


334.5 






16 


91 


(0 


07) 


14 


72 


(0.03) 


RCC 


20041130 


340.4 


19 


04 (0.06) 


17 


16 


(0 


04) 


14 


87 


(0.04) 


RCC 


20041211 


350.5 


19 


56 (0.08) 


17 


60 


(0 


05) 


15 


23 


(0.06) 


RCC 


2004 12 22 


361.5 


18 


87 (0.11) 


17 


24 


(0 


08) 


14 


91 


(0.05) 


RCC 


2004 12 30 


370.3 


18 


92 (0.05) 


17 


24 


(0 


06) 


15 


01 


(0.05) 


RCC 


2005 01 06 


377.3 






16 


95 


(0 


06) 


14 


79 


(0.07) 


RCC 


2005 01 10 


381.3 


18 


74 (0.06) 


17 


08 


(0 


08) 


14 


91 


(0.10) 


RCC 


2005 0111 


382.3 






17 


04 


(0 


07) 


14 


84 


(0.07) 


RCC 


2005 01 20 


391.3 






17 


35 


(0 


06) 








IAC-80 


2005 02 01 


403.3 






17 


67 


(0 


04) 


15 


34 


(0.10) 


RCC 


2005 02 05 


407.3 






17 


43 


(0 


09) 


15 


23 


(0.05) 


RCC 


2005 02 07 


409.3 


18 


93 (0.09) 


17 


48 


(0 


07) 


15 


09 


(0.05) 


RCC 


2005 03 21 


451.3 


19 


06 (0.07) 


17 


18 


(0 


03) 








IAC-80 


2005 10 04 


648.5 


20 


35 (0.10) 


18 


65 


(0 


06) 


16 


00 


(0.08) 


RCC 


20051005 


649.5 


20 


43 (0.06) 


18 


58 


(0 


04) 


16 


14 


(0.05) 


RCC 


20051009 


653.6 














16 


39 


(0.06) 


RCC 


2005 10 10 


654.6 














16 


16 


(0.04) 


RCC 
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Table 1- 


— Continued 






Date 




V 




Ic 


Telescope 


2005 10 15 


659.7 






16.03 (0.07) 


RCC 


2005 10 19 


663.6 






16.29 (0.05) 


RCC 


2005 10 28 


672.6 


21.34(0.16) 


19.49 (0.08) 


17.02 (0.08) 


RCC 


2005 10 31 


674.6 


21.55 (0.16) 


19.66 (0.05) 


17.17(0.03) 


RCC 


20051101 


675.6 


21.75 (0.08) 


19.83 (0.04) 


17.30 (0.04) 


RCC 


20051117 


692.5 






17.55 (0.08) 


RCC 


20051119 


694.6 






17.70 (0.10) 


RCC 


2005 12 08 


713.4 






17.88 (0.10) 


RCC 


2005 12 10 


715.4 






17.95 (0.10) 


RCC 


2005 12 11 


716.4 






17.63 (0.10) 


RCC 


2005 12 13 


718.5 






17.84 (0.10) 


RCC 


2006 02 01 


768.3 






19.49 (0.10) 


RCC 



^2,453,000+ 
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2.2. Near-infrared imaging, photometry and polarimetry 

Near-infrared J, H and band observations were carried out at 10 epochs in the period 2004 
March - 2006 March. The data were obtained using two instruments: LIRIS, instahed on the 4.2 m 
Wilham Herschel Telescope (WHT) at the Observatorio del Roque de Los Muchachos, and CAIN-2, 
installed on the 1.52 m Carlos Sanchez Telescope (CST) at the Teide Observatory. Table [2] lists the 
epochs of all near-infrared observations. In addition, polarization images through a J filter were 
obtained during the night of 2004 October 28, using the instrument LIRIS. 

LIRIS is a near-infrared camera/spectrograph attach ed to the Cassegrain focus of the WHT 



(jAcosta-Pulido et al.l l2003l : iManchado-Torres et al.l 120031 ) . The data were obtained during the 



commissioning period (2004 March) plus some nights corresponding to the guaranteed time of the 
instrument team. LIRIS is based on a Hawaii-I detector with an image scale of 0'.'25 per pixel and 
4'.2 X 4'.2 field of view. Observations were always performed by taking several dithered exposures 
around the position of V1647 Ori in order to ensure proper sky image subtraction. For each 
dither point we took several frames using short exposure time to avoid saturation, which were later 
averaged to obtain the final reduced image. The minimum exposure time was 1 s for individual 
exposures, the total integration time per filter was about 5 min. During the first observing period 
the individual exposure times were optimized to detect the nebula, which nearly saturated the 
detector at the position of the illuminating star. 

The infrared camera CAIN-2 at the 1.52 m TCS is equipped with a 256x256 Nicmos 3 detector, 
which provides a pixel projection of 1" with its wide optics configuration. The same dithering 
technique was used as in the case of LIRIS observations. Total integration time was usually 10 min 
in all filters (J, H and K^), splitted in many exposures of 0.5-3 s, the shortest corresponding to 
filter Ks and the longest to filter J. 



2.2.1. Data reduction and photometric calibration 

Data reduction process involved the following steps: sky subtraction, flat-fielding and coaddi- 
tion of all frames taken with the same filter. The sky image was obtained as the median combination 
of all frames, masking regions occupied by bright sources. The final image was produced using the 
standard "shift-and-add" technique including rejection of outlier pixels. The LIRIS images were 
reduced using the package 'liris_qr developed within the IRAF environment. The CAIN-2 data 
were also treated with IRAF tasks. The J and Kg images obtained on 2004 November 4 can be 
seen in Fig. [3J 

The instrumental magnitudes of the central star were extracted using aperture photometry 
since the contribution of the nebula in the near-infrared becomes negligible compared to it. Based on 
the measured radial profiles we estimate that the nebula contribute to the central star photometry 
at most 15% in the J band, but less than 5% in the Ks band. Due to the high brightness of 
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.-0 



30 20 10 -10 -20 

OffRA (arcsec) 



Fig. 3. — Near-infrared images of McNeil's Nebula obtained on 2004 November 4 with LIRIS 
on WHT. North is at the top and east to the left. The color scaled image represents the K<j- 
band emission and the contours represent the J-band emission. The inset is a zoom on the star 
southwest of V1647 Ori, demonstrating that the star is a double system. The o pen diamond at 
(Ag, A .6) ~ (16", 38") marks the position of the bright infrared source detected by 
(|200J) in Spitzer images. 



Muzerolle et al 



- 12 - 



V1647 Ori some LIRIS images entered the non-linear regime at the position of the star. In these 
cases we estimated the flux of the star by fitting a PSF model to the non-saturated parts of the 
stellar profile using the daophot package in IRAF. 



For the photometric calibration we used the 2MASS catalogue (ICutri et alJ 120031 ). In general, 
five or six 2MASS stars were found within the sky region covered by our images. We have determined 
the offset between the instrumental and the calibrated 2MASS magnitudes by averaging al l available 
stars, after removal of deviant sources (identified as likely variable stars bv lSemkov II2006I ). Results 
of the near-infrared photometry for V1647 Ori are listed in Tab. [21 and the light curves are shown 
in Fig. [2l The standard deviation is used as an estimate of the measurement error and reported 
in Tab. [2j We cannot determine a reliable color term correction using our reduced number of 
reference stars. The 2MASS star J054611162-0006279 was always included as a comparison star, 
since it shows calibration offsets clos e to the average in all occasions. We note that this source 
was the only comparison star used by IWalter et al.l (120041 ) for the photometric calibration of their 
near infrared data. The colors of this star are very similar to those of V1647 Ori which reduces the 
importance of color terms when determining the photometrical calibration. 



2. 2. 2. Polarimetry 



The polarization images were obtained thanks to a double Wollaston prism (WeDoWo) (jOliva 



19971 ). located in the LIRIS grism wheel. This device produces four different images on the detector 
from a rectangular field of 4 x 1 arcmin^. Each image corresponds to a different linear polariza- 
tion angle, namely 0, 90, 45 and 135 degrees. After distorsion correction, proper alignment and 
combination of these images, the degree and angle of linear polarization can be determined. For 
V1647 Ori and its nebula we obtained an exposure of about 1800 s, divided in individual frames of 
15 s. A small dither pattern of 3 points was performed around the source, interleaved with mea- 
surements at an offset position in order to determine a proper sky emission, given the extension of 
the source. The data reduction process is very similar to that followed in normal imaging mode, the 
'liris_qr package was used. Images of the linear polarization degree will be presented and discussed 
in Section [5T3l 



2.3. Near-infrared spectroscopy 

Near-infrared spectrograms of V1647 Ori were obtained using LIRIS at five different epochs, 
four of them during the first year of the outburst. Five spectra cover the wavelength interval of 
0.9-1.4^m {ZJ bands) and three others cover the 1.4-2.4^m (HK) wavelength band. The log of 
the observations is presented in Tab. El 

Observations were performed following an ABBA telescope nodding pattern both in the ZJ 
and HK bands. In order to reduce the readout noise, the measurements were done using multiple 
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Table 2. Near infrared photometry of V1647 Ori 





J JJ 


J 


tl 




Telescope 


2004 03 04 


069.4 


11.13 (0.10) 


9.09 (0.07) 


7.48 (0.09) 


LIRIS 


2004 03 21 


086.3 


11.00 (0.12) 


8.70 (0.19) 




TCS 


2004 09 01 


249.6 


11.12 (0.10) 


9.17(0.03) 


7.66 (0.01) 


TCS 


2004 10 05 


283.7 


11.06 (0.07) 


9.07 (0.14) 


7.76 (0.05) 


TCS 


2004 10 17 


295.7 


11.10 (0.09) 


9.09 (0.14) 


7.62 (0.04) 


TCS 


200410 21 


299.8 


11.17(0.08) 


9.19(0.05) 


7.84 (0.08) 


TCS 


20041104 


313.7 


11.20 (0.06) 




7.77(0.06) 


LIRIS 


2006 01 14 


749.6 


14.27(0.09) 


11.48 (0.08) 


9.41 (0.19) 


LIRIS 


2006 02 06 


770.8 






9.78 (0.02) 


TCS 


2006 03 31 


826.5 


14.19(0.13) 


11.55 (0.05) 


9.57 (0.08) 


TCS 


2006 05 05 


861.4 


13.92 (0.08) 






LIRIS 


2006 09 09 


987.8 




11.59 (0.05) 


9.82 (0.06) 


LIRIS 



2,453,000+ 



Tabic 3. .JoTirnal of Ncar-infrarcd spectroscopy 



Date 




Slit 


Band 


^exp 


Band 


^exp 


Airmass 










(s) 




(S) 




2004 03 08 


073.4 


1 


ZJ 


400 


HK 


12 


1.3 


2004 11 04 


313.7 


0'.'75 


ZJ 


800 


HK 


432 


1.2 


2005 01 24 


395.4 


0^'75 


ZJ 


200 






1.2 


2005 03 25 


455.4 


1 


ZJ 


100 






2.0 


2006 05 05 


861.4 


1 


ZJ 


1200 


HK 


400 


4.1 


2006 09 09 


987.7 


a'75 


ZJ 


1200 


HK 


900 


1.4 



2,453,000+ 
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correlated readout mode, with 4 readouts before and after the integration. We used a sht width 
of CK'TS or 1", depending on the seeing conditions, which yielded a spectral resolution in the range 
R = 500-660 and 550-700 in the ZJ and HK spectra, respectively. The wavelength calibration was 
provided by observations of Argon and Xenon lamps available in the calibration unit at the A&G 
box of the telescope. In order to obtain the telluric correction and the flux calibration, nearby AOV 
or G2V stars were observed with the same configuration as the object. 

The data were reduced and calibrated using the package 'liris_qr. Consecutive pairs of AB 
two-dimensional images were subtracted to remove the sky background, then the resulting images 
were wavelength calibrated and flat-fielded before registering and coadding all frames to provide 
the final combined spectrum. One dimensional spectra were extracted with the IRAF 'apall' task. 
The extracted spectra were divided by a composite spectrum to eliminate telluric contamination. 
This composite spectrum was generated from the observed spectra of the calibration star, divided 
by a stell ar model and conv olved to our spectral resolution. We used the publicly available code 



'Telluric' (jVacca et al.ll2003l ) to perform the correction when AO stars were observed. Differences in 
the strength of telluric features likely due to mismatch of air masses and variation of atmospheric 
conditions between observations of the object and the reference star were taken into account using 
Beer's law. The IRAF task 'telluric' was used in this step. The fiux calibration was carried out 
normalizing to the J, H, and magnitudes obtained in our near-IR photometry. We show in 
Fig. m the flux calibrated infrared spectra for those observing dates when the full range from 0.9 to 
2.4/xm was observed. 



2.4. Optical spectroscopy 

Medium-resolution optical spectra were taken on three nights: 2004 October 28, November 1, 
and December 24, with the Double Beam Spectrograph mounted on the 2.3 m telescope of the Siding 
Spring Observatory (Australia). We used 1200 mm~^ gratings in both arms of the spectrograph, 
while the exposure time was set to 20 minutes. The slit width was 2" and the instrument's rotator 
kept it in the parallactic angle. The image scale was O'.'O/pixel in the direction of the dispersion 
and 3'.'6/pixel perpendicular to it. The observed spectral ranges and resolutions were as follows: 
blue - 4800-5240 A, A/AA=8300; red - 5820-6780 A, A/AA=6500. Ah spectra were reduced with 
standard IRAF tasks, including bias and flat field corrections, aperture extraction and wavelength 
calibration. 

In the blue spectral range w e observed a barely visible continuum, and could not detect the 



absorption spectrum observed by iBriceho et al.l (|2004l ) in 2004 February. These blue spectra were 



excluded from the further analysis. In the red range we observed the Ha emission line on a fea- 
tureless continuum. The equivalent widths of the stellar Ha line determined from our observations 
are 32, 34 and 27 A for the dates 2004 Oct 28, Nov 1 and Dec 24, respectively. 
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Fig. 4. — Flux calibrated near-infrared spectra of V1647 Ori obtained with LIRIS on WHT. 
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3. Brightness variations during the V1647 Ori outburst 
3.1. Long term evolution 



In order to examine th e flux evolution duri ng the whole outburst history we combined our Iq 
measurements with those of lBriceho et al.l (|2004l ). which covered the beginning of the eruption. The 
complete /c-band light curve is displayed in Fig. [5l Following an initial brightening of ~ 4.5 mag 
between 2003 November and 2004 February (jBricefio et al.ll2004l ). the star faded only moderately 
during the next one and a half years. This gradual decrease of brightness was observable in all 
photometric bands with similar rates (Fig. [2]). Then, between 2005 October and November the 



brigh tness of the star suddenly dropped b y more th a n 1 rn ag at all optical wavelen gths (jKospal et al 



20051 ). This decay was also monitored bv ISemkov I (|2006l ) and lOjha et al.l (|2006l ). The rapid fading 
continued in the next months, although the rate was slightly slower after 2005 mid-November. By 
2006 February, some 800 days after the onset of the eruption, V1647 Ori was very close to the 
reported pre-outburst state. 

The complete light curve of the outburst can naturally be divided into three phases: the initial 
rise, the "plateau", and the final decay. In order to quantify the rate of brightness evolution, 
we fitted straight lines to the data points in Fig. [5] for each phase separately. The results are 
— 1.3 mag/month, 0.04 mag/month and 0.8 mag/month, respectively. It is remarkable that the 
pace of brightness change in the rising and decaying phases is similar. 

The fact that the fading rate of V1647 Ori is apparently wavelength independent is different 
from the case of another eruptive star V1057 Cyg, whose decline rate decreas ed monotonically from 
the u ltraviolet to the mid-infrared during the years following the outburst (jKenyon Hartmann 
). However the decline of V1057 Cyg was most pronounced in the B band, which was not 



1991 



observable in the case of V1647 Ori. 



3.2. A periodicity of 56 days 



The light curves in Figs. [2] and [5] show that, in addition to the slow decline during the high 
plateau, short time scale variations can be seen in each photometric band. Definite brightness min- 
ima can be recognized in each optical band e. g. on 2004 December 10 (JD 2,453,350) and on 2005 
February 1 (JD 2,453,403). Similar depressions can be seen around 2004 March 7 [JD 2 ,453,072) 
and 2004 April 12 (JD 2,453,108) in the light curves presented by Iwalter et al.l (j200J). Visual 
inspection suggests a possible periodicity in the light curve. In order to study this phenomenon 
quantitatively , a period sear ch analysis was conducted, using the programme package MUFRAN 
developed bv iKollathl (|l99d ). The MUFRAN software is based on the Fou r ier tr ansform. For a 
recent review of the time-series analysis of variable star data, see iTempleton I (j2004l ) and references 
therein. 



Our Rq and Ic data listed in Tab. [T] were supplemented with data taken from lSemkov (2004, 
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20061 ) in the Rc and Ic filters. The imported data were zero offset by 0.42 (Re) and 0.17 magnitudes 



(/c), in order to reach a common average magnitude scale with our data. The most likely reason 
of such difference is the use of different methods to extract the stellar flux. Only the time interval 
between 2004 Aug 18 and 2005 February 11 (JD 2,453,235-JD 2,453,413) was considered because 
outside this interval the longer time-scale variation may falsify the result of the period analysis. 
Due to the evidence of a long term brightness variation, allowance was made for the monotonous 
dimming: a straight line was fitted to each dataset and subtracted to form magnitude differences 
with respect to this line. The period search was performed separately for data in each photometric 
band. In the case of finite dataset and unevenly sampled data some false frequencies (aliasing) 
can appear in the Fourier power spectrum. These alias frequencies are centered on the real signals 
offset from those peaks as indicated by the spectral window function of the Fourier transform. The 
temporal distribution of the photometric data of V1647 Ori is favourable and the spectral window 
(see the left panel of Fig. [6|) indicates a negligible aliasing caused by the data sampling. Here 
we only show the spectral window function corresponding to the Iq filter dataset, and note that 
the shape and the alias structure is practically the same in the case of the spectral window for 
the Rc dataset. The power spectra typically contain a number of peaks (see middle panel of Fig. 
E]). We show in this figure only the lowest frequency part of the power spectrum, the rest of the 
power spectrum only contains much lower peaks than seen here. The highest peak appears at the 
frequency of / = 0.0179 cycle/day (equivalent to a period of 55.87 days) in the Iq data. The S/N 
ratio is slightly above 3 in the region of the / = 0.0179 c/d peak, which gives significance to this 
peak. The S/N value was determined as the ratio of the amplitude of the main peak and the average 
amplitude present in the power spectrum after prewhitening with this frequency. Moreover, a peak 
can be considered as physically meaningful, if it appears in the power spectra of each photometric 
observational series. The corresponding peak is situated at the frequency of 0.0163 c/d in the 
analysis of the shorter series of the Rq filter data. The phase curve of the Iq residual photometric 
data folded on the 55.87 day period is shown in the right panel of Fig. [6l The Rc band folded light 
curve is qualitatively similar to that shown in Fig. [6l It is worth mentioning that the negligible 
phase difference between the Rc and Ic folded light curves also implies a common origin of this 
periodicity. This means that the ~ 56 day periodicity is confirmed. We have plotted in Fig. [7] the 
portion of the Ic light curve which was fitted with the sinusoidal curve plus a monotonic decrease. 
The time interval covered by our data is 200 days, which is only about 4 cycles of the possible period 
of 56 days. The periodic behaviour is clearly seen in the curve, although a sinusoidal variation is 
regarded only as a first and crude approximation. 



3.3. Hour scale variations 



Ojha et al.l (120051 ) claimed weekly variations of the optical fiux of V1647 Ori. iGrosso et al 



reported on a periodic variation in the X-ray fiux of 0.72 days, about 17 hours. We have 



also searched for very short term variability in our optical data. We compared 5 consecutive Ic- 
band observations of V1647 Ori obtained on the same night (2004 November 30) during a period 
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Fig. 5. — /c-band light curve of V1647 0ri covering the complete 2-v ears outburst period. Filled 
stars: our observa tions reported in Tab. [H open circles: data from iBricefio et al.l (|2004l ): open 
squares: data from lOjha et al.l (j2006l ). The uncertainties are comparable with symbol sizes. 
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Fig. 6. — Results of the period search. Left: The spectral window function for the Ic photometric 
sample; Middle: A part of the Fourier power spectrum of the Ic data. The first high peak at the 
frequency of 0.0179 c/d corresponds to a period of 55.87 days. Right: The Ic residual phase curve 
folded on the period of 55.87 days. 
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of 2.5 hours with the 1 m RCC telescope. The first four measurements are identical within their 
formal photometric uncertainties (~ 0.03 mag), but the fifth measurement deviates about 0.2 mag 
from the others. This figure may be the characteristic variation expected on hour scale. 



3.4. The possible origin of variability 



In order to study further the mechanisms causing the initial brightening and the final fading, 
we plotted in Fig. [S] optical and near-IR color-magnitude diagrams, and looked for relationships 
between bright ness and color ind ex variations. We also marked in Fig. [5] the standard interstellar 
reddening law ( Cohen et al.lll98l1 ). 

Reipurth fc Aspin I ^20oi ) noted that the star mo ved along the exti nction path in the J-H 



vs. H-K diagram during the initial brightening phase. lOjha et al.l (120061 ) demonstrated that the 
final fading of the object in November 2005 also followed the same path in the opposite direction. 
However the variation of the near-infrared and visible colors is not compatible with simple dust 
obscuration. It can be seen from Fig. [8] that the amount of extinction needed to explain the color 
variation in the near-infrared is much larger than the one required in the visible. In order to explain 
the amplitude of the outb urst a color independe nt flux increase has also to be assumed. A similar 
conclusion was reached bv lMcGehee et al.l (j2004l ). 



In order to explain the flux variations during the plateau phase we looked at any possible 
relationship between brightness and color index variations (Fig. E ). A loose trend that the star 
tends to be redder when dimmer can be observed in the diagrams. lEiroa et al.l (l2002l ) successfully 
explained the light variations of several young stars by obscuration by dense slabs within the 
circumstellar disk. Thus obscuration by circumstellar dust structures orbiting at a distance of i? = 
0.28{M/Mq)'^/^AU, implied by our ~ 56 day period, may explain at least part of the photometric 
variations during the plateau phase. 




200 



250 



300 350 
JD (2,453,000 + ) 



400 



Fig. 7. — The Ic light curve with the fitted sinusoidal curve overimposed. The dimming trend was 
also included. The quasi-periodic variation of the light curve is clearly detected. 
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3.5. Inclination of the star nebula axis derived from the brightness evolution 

McNeil's Nebula shines due to the scattered ligh t of V1647 Qri, t heref ore its observed properties 



are closely related to the brightness of the star. iBricefio et alj (|2004l ) examined the temporal 
evolution of two well-defined positions of the nebula (marked as B and C in their fig. 2), and 
concluded that their light curves followed the brightening of the star with delays of about 14 days 
(Pos. B) and of more than 50 days (Pos. C). The angular distance of 36" between the star and 
Pos. C implies a delay of 85 days. Since these delays carry information on the incl ination of the 



nebu la (see Fig.[9j), we reanalysed the three /c-band light curves (from Fig. 3 in IBricefio et al 



20041 ). supplemented with our measurements obtained in the period 2004 February- April, in order 
to find the accurate time shifts between them. First we interpolated the light curves for a time 
sampling of 0.25 days. The resulting light curves of V1647 0ri, Pos. B, and Pos. C are called A{ti), 
B{ti) and C{ti), respectively, where ti is the time [i = \...N). Then we calculated a function called 
D{s) as follows: 

N 



D{s) = Y,A{ti)-C{ti-s) 



The quantity s which results in the lowest value of D{s) gives the time shift between the light 
curves of A and C. The shifts between A and B can be calculated similarly. In this manner we 
obtained a time shift 12.5ib2 days between V1647 0ri and B and 50ib3 days between V1647 0ri 
and Pos. C . 

The measured time delay between the light curves of V1647 Ori and Pos. C, together with 
the projected separation of the two objects results in an angle of 6II2 degrees between the axis of 
the nebula and the line of sight. Because the clump lies inside the polar stellar wind cavity of the 
star, this angle is a good estimate for the inclination of the rotation axis of the star (see Fig. [9]). 
We estimated the uncertai nty of the inclination der ived this way as follows. We assumed that the 



Herbig-Haro object HH23 (jEisloffel &: Mundtlll99?1 ) lies exactly along the rotation axis of the star 



Then we assumed that the projected angle between the star-Pos. C and star~HH 23 directions, 
11°, is the same as the difference in their tilts to the line of sight. Thus the total uncertainty of 
the inclination derived from the time delay of the brightness variation of clump C is ±14°. 

The time delay between the light curve features of V1647 0ri and Pos. B results in 65°ib8°. 
The inclinations derived from positions B and C are compatible within the errors. Due to the larger 
opening angle of the conical nebula at small distances from the star, position B is less suitable for 
estimating the inclination of the stellar rotation axis. 



This result apparently conflicts with lRettig et al. I (|2005l ) finding that the disk of V1647 Ori is 



seen nearly face-on, because the absorption line of the cold CO due to the outer, flared regions of 
the disk cannot be observed in the high resolution infrared spectrum of the star. But if the thickness 
H of the disk at a distance R from the star does not exceed H/R ~ 0.50, then the absorption lines 
produced by the outer disk is not expected to be observed in the spectrum of the star. 
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4. Spectral evolution 



Our near-infrared spectra, normalized to the continuum, are displayed in Fig. [10] (ZJ band) 
and Fig. [11] (HK band). The most prominent features are the Pa/3, Pa7, and Br7 lines of HI in 
emission, the Hel 1.083 ^m line most of the time seen in absorption, the Mgl 1.50 fj,m line, and the 
rovibrational Aiy = +2 transitions of CO at 2.3-2.4 fim in emission. These spectral features are 
present in all spectra, though significant changes can be recognised in the line strengths during our 
observational period of 30 months . In th e following w e analy se th ese spectral variat ions, including 
also data sets from IWalter et al.l (|200J) , iGibb et al.l (|200a) and lOjha et al.l (|200a) , which cover 
shorter periods. We note that our last spectra were taken in 2006 May and Sep when the outburst 
was over, thus probably they are representative of the quiescent phase of the object. 

Table H] lists the equivalent widths and fluxes of the infrared lines. The near-infrared line 
fluxes were computed from the absolutely calibrated spectra, with a typical estimated uncertainty 
of about 15% in the brightest lines to 25% in the faintest ones. Figure [12] shows the time evolu- 
tion of the line fluxes. In addition, we plotted the Ha and Call A 8542 A l i nes, w here line fluxes 



were d erived by combining the equivalent widths published by 1 Walter et al.l ([200J) and lOjha et al 



(|2006l ). and Ic and Rc magnitudes taken from Tab. [T] or from the literature for the same day as 
the spectroscopic data. The typical uncertainty of these latter flux values, derived from the errors 
of the magnitudes and equivalent widths, is about 6 percent. We noticed inconsistencies when 
comparing our line fluxes with those published by other authors for close observing dates. For this 
reason we recomputed their line fluxes combining the reported EW with photometric data. These 
are the values plotted in Fig. [12] 



Hydrogen lines. Figure [12] shows that the line fluxes of Pa/3, Br7 and Ha declined slowly during 
the period of observations. The Pa7 shows a similar trend, but with higher scatter. 



In many young stellar objects the origin of HI lines is related to accretion processes. iMuzerolle et al 



(jl998al ) established empirical relationships between the accretion rate and the luminosity of Pa/3 
and Br7 lines. The data in Fig. [12] clearly indicate parallel decrease of the strongest hydrogen emis- 
sion lines during the whole outburst, which wo uld suggest a decreasing accretion rate. Application 
of the relationships of iMuzerolle et al.l (|l998al ) with the assumptions that V1647 Ori is located at 
400 pc from the Sun; suffers an extinction of Ay = 11 mag; has a mass of 0.5 M© and an accretion 
radius of SR©, results in an accretion rate about 5 x 10~^MQ/yr in 2004 March, and an order of 
magnitude lower value in 2006 May. 

The Paschen lines of HI showed P Cygni profile with a blueshifted absorption o i ily in 2004 
March, whereas the absorption part disappeared in later observations. iGibb et al.l (|2006l ) also 
reported the same vari ations in the profiles of the HI lines. Similar changes of the Ha and Call 
lines were reported by lOjha et al.l (|2006l ). These results indicate that at the beginning of the 
outburst the source of the HI lines was at least partly stellar wind, whose strength decreased 
during the outburst. As a consequence, the drop of accretion rate was probably smaller than we 
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Fig. 8. — Relationship between the color indices and magnitudes of V1647 Ori during the monitoring 
period. Left: V vs. V — Iq] Right: J vs. J — Kf^ . The filled symbols represent our measurements, 
and the open symbols are take n fromlOiha et alj (|2006l ). The dotted lines indicate the slope of the 
normal interstellar reddening (ICohen et al.lll98lh . Solid lines connect consecutive observed points. 



Table 4. Equivalent width and flux of near-infrared lines. 



Line 


A (/im) 


2004 Mar 


2004 Nov 


2005 Jan 


2005 Mar 


2006 May 


2006 Sep 






EW^ 


Flux^ 


EW 


Flux 


EW 


Flux 


EW 


Flux 


EW 


Flux 


EW 


Flux 


Pa5 


1.005 


-0.8 


1.2 


-2.5 


1.2 


-2.1 


0.5 


-1.7 


0.4 






-12.4 


0.26 


Hel 


1.083 


9.18 


-5.9 


4.3 


-3.0 


3.3 


-1.1 


5.2 


-1.8 


-11.6 


0.3 


-17.0 


0.56 


Pa7 


1.094 


-3.0 


1.9 


-2.4 


1.7 


-2.7 


0.9 


-2.2 


0.8 


-14.8 


0.5 


-13.4 


0.50 


Pa/? 


1.282 


-10.7 


14.4 


-8.5 


10.1 


-12.5 


8.1 


-15.1 


8.2 


-14.9 


1.4 


-16.7 


1.51 


Mgl 


1.504 


-3.7 


7.6 


-2.8 


4.4 


-5.9 


6.6 


-6.7 


5.2 






-1.3 


0.23 


Bi5 


1.945 


-2.9 


10.5 






















Hel 


2.058 


< 0.8 


< 3.0 


0.8 


-2.5 














-1.7 


0.77 


Br7 


2.166 


-5.2 


22.0 


-4.8 


15.1 










-6.9 


4.0 


-6.7 


3.31 



^Units A 

^Units lO-^'^Wm-^ 




Fig. 9. — Geometry of McNeil's Nebula, deduced from the time delay between the light curves of 
V1647 Ori and Clump C. 
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Fig. 10. — Near-infrared {ZJ band) spectra of V1647 Ori obtained with LIRIS, normalized to the 
continuum. 
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derived in the previous paragraph. 



Walter et alj (j2004l ) and lOiha et al.l (j2006l ) suggested a possible periodic behaviour in the Ha 



equivalent width with a period of roughly 50 days. Since this period is similar to the one found 
in our Iq and Rc light curves (55.87 days, see Sect. 3.2), it is tempting to speculate about a 
connection between the two cyclic phenomena. In order to check it, we folded all available Ha 
equivalent width data with the period of the optical light curve, and constructed a phase curve. It 
was found that the EW(Ha) and Ic values oscillate in opposite phase, i.e. the equivalent width of 
the Ha line is the largest when the central object is the dimmest. Since the equivalent width is the 
ratio of the true line flux to the continuum level, it is very likely that the periodic variations of the 
Ha EW are not intrinsic but due to the periodically changing continuum flux. 



Helium lines. The evolution of the Hel 1.083 /im line tells an interesting story. It was seen 
in blueshifted ab sorption during t he whole outburst, although its strength decreased with time 
(Fig. [T2| see also iGibb et al.l 120061 ). The flux decreases quickly towards the end of 2004, and later 
decreases again in 2006 May. In our spectra of 2006 May and Sep, taken in the quiescent phase 
this line was detected in emission, without any apparent velocity shift (Fig. [TO]) . The transition 
of the Hel 1.083 ^m line from absorption to emission is a remarkable event, not seen in any other 
line in our spectra. This line also shows a blueshifted (~ 500km s~^) absorption component in the 
spectrum taken in 2006 Sep. 

In young stellar objects the Hel 1.083 /xm line is usually attributed to hot winds. The profile of 
the Hel absorption line o bserved in V1647 Or i was similar to those seen in a large sample of strongly 
accreting T Tauri stars (lEdwards et al.ll2006l ). where the wh ole velocity raiige of t he Hel wind can 
be observed against the stellar continuum. In these cases lEdwards et al.l ([200a) argue that the 
wind traced by the Hel 1.083 //m line originates close to the star and not from the inner part of the 
circumstellar disk. This may also be the case for the V1647 Ori system. The hard X-ray emission 



from V1647 0ri, observe d by the X-ray observatories Chandra and XMM-Newton ([Kastner et al 



20041 : iGrosso et al.ll2005l ). indicates the presence of a magnetic reconnection ring located at 1-1.5 
stellar radii above the stellar surface in the plane of the disk. We speculate on that the source 
of the wind is the magnetic reconnection ring. This hy pothesis is supported by the fact th at the 
strongest Hel hne, observed on 2004 March 8, (Fig. [T2l Ivacca et al.ll2004l: bibb et ahl bood). was 



preceded by a high X-ray emission event on 2004 March 7, (|Kastner et al.ll20o'i K indicating a direct 



connection between the two phenomena. The drop of the flux of the Hel 1.083 /um line during the 
outburst suggests weakening of the wind. 



Edwards et al.l (j2006l ) proposed that there may be two types of the hot wind, producing ab- 



sorption or emission of the Hel 1.083 /im line. The transition of this line from absorption to emission 
at the end of the outburst indicates that the nature of the wind in V1647 Ori changed between the 
two types of hot wind. 



We also detected the Hel 2.058 /im line in absorption both in 2004 November and 2006 May. 
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The apparent increase of its strength in the latter spectrum (taken during quiescence) is due to the 
lower continuum. 



Other spectral features. In Fig. [12] we also plotted the flux of the Mgl 1.505 //m emission line. 
The feature was detectable during the whole outburst, showing hints for a slowly decreasing trend 
in flux. In the spectrum of May 2006, during the quiescent phase, the line was not visible; the 
derived upper limit in Tab. H] indicates that the Mgl 1.505 /im emission is weak in the quiescent 
phase. 

The last panel of Fig. [12] shows the flux variation of the Call 8542 A emission line based on 
literature data. Th is middle component o f the Call triplet was found to be an excellent tracer of 
the accretion rate (jMuzerolle et al.lll998bl ). It can be seen that, contrary to all the other lines in 
Fig. [T2] the flux of the Call 8542 A line was constant during the period of the observations (there is 
no available measurement from the quiescent phase). A reason for the discrepancy may be that the 
hydrogen and cal cium lines originate from regions where the physical conditions are very different. 
Ojha et al.l (l2006l'l noted that the den sity of the Call region is very high and very optically thick. 



According to iMuzerolle et al. I (120011 ) optically thick lines are less reliable tracers of the variable 
accretion. The accretion rate, derived from its luminosity, is about 2 x W^^Mq j yr, close to the 
values obtained from the hydrogen lines for the high state of the outburst. 

The CO band head features a t 2.3-2.4 jxm. were observ e d in strong emission during the whole 
outburst period (Fig. [TTl see also iReipurth fc Aspin II2004I : IVacca et al.ll2004l : iRettig et al. II2005I : 



Gibb et al.ll2006l ). However, the CO feature is not present in our quiescent spectrum suggesting 



that the appearance of this band is related to the outburst. 



5. The nebula around V1647 Ori 



McNeil's nebula is a reflection nebulosity which scatters the light of its illuminating star, 
V1647 Ori. The optical and near-infrared morphology of the nebula can be seen in Fig.[T] and 
Fig. [31 respectively. Their study provides information not only on the neighbouring interstellar 
medium but als o on the outbur s ting star. In the follow ing we describe the morphology of the 
nebula (see also lOjha et al.ll2005l : iReipurth fc Aspin II2004I ). 



Comparing the optical and near-infrared images one immediately notices that the shape of the 
nebula looks rather different at optical (Fig. [1]) and at infrared wavelengths (Fig. [3]). In the optical 
images the nebula shows a conical shape, extending up to a distance of ~ 50" (~ 22000 AU at the 
distance of the nebula) to the north. The diffuse emission area appears to be delimited by two 
bright regi ons indicating an op ening angle of ~ 50°. The north-west side is marked by region B 
(named by lBriceho et al.ll2004l ). a plume located at a distance between 5 and 10". The no rth-east 
rim appears more extended and is delimited by region C (also named by lBriceho et al.ll2004l ). which 
is located around 35" from the star and is elongated perpendicularly to the radial axis from the 
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star. No southern extension of the nebula can be recognized in the optical images. 

In contrast to the optical, the infrared morphology indicates a much more compact structure, 
though the cometary shape can still be seen with a sharp curved tail along the north-east side. 
The nebula appears to be more extended in the J band, exhibiting diffuse emission both at the 
north-east and west sides. It appears roughly spherical in the band, although the structure is 
flattened perpendicularly to the nebula axis; the largest size of the structure is ~ 20". This size 
corresponds to a diamet er of ^ 9000 AU a ,t the distance of Orion B cloud, which is slightly smaller 
than the size derived bv I Johnstone et al.l (|200ll ) for the matching submillimeter source OriB N 55 , 
about 13500 AU. The southward extension is also clearly seen in the band. 



5.1. Color variations across the nebula 

In Fig. [13] we show the resulting optical (V— R) and near-infrared (H— K) maps. In addition 
we obtained cross-cuts along the N-S direction passing through the star (Fig. [T4l) . 

Some distinguishable regions can be recognised in the V — Rc map: i - the bluest colors of 
the nebula are registered towards the North end of the nebula, its color is around V — Rc ~ 
and Rc — Ic ~ 0.3 at 60" from the star. Assuming that the light path along this direction is 
relatively cleaned of dust, we may observe the colors of the closest vicinity of the star, though 
somewhat modified by scattering and foreground extinction. These colors are very close to those 
of an intermediate A type star, indicating the pres ence of a hot central area probably related to 



the outburst of this otherwise cold, late type star (jMcGehee et al.ll2004l ). ii - the plume located 
at a distance between 5 and 10" to the north-west (position B as defined by lBriceho et al.ll2004l ) 
appears redder than the rest of the nebula, in - a compact region of strong reddening is observed 
towards the direction of V1647 Ori. 

In the H — Kg map the reddest part of the nebula can be seen to the south of the illuminating 
star, outlining a flattened D-shape structure, whose elongation is nearly perpendicular to the axis of 
the optical nebula. The nebula appears very red along the E-W direction across the star (J—H ~ 2 
and H — Kg ~ 1). When crossing the central part towards North, the nebula changes rapidly to 
bluer colors {J — H and H — Kg ~ 0.5), although the extension of the nebula is smaller than in 
the visible range. The infrared colors are redder at the South compared to the North values at the 
same radial distance, indicating a larger extinction along the line of sight in the South part, or part 
of emission from the dust envelope. 

The observed color distributions of McNeil's nebula can be interpreted in terms of a thick 
circumstellar envelope which intersects the light path from the central parts of the system to the 
scattering grains in the outer envelope. The envelope includes conical outflow cavities perpendicular 
to the disk plane, whose walls scatter the starlight. These cavities had been cleaned up by previous 
eruptions, and allow the light to escape and reach the northernmost regions of the nebula. The 
North-South asymmetry is likely caused by the inclination and the different amount of extinction 
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along our line of sight. A similar geor netry is assunied in theoretical models of infrared re flection 
nebulae around young stellar objects (jLazareff et al.lll99d : iFischer. Henning &: Yorkdll996l ). 



5.2. Extended line emission 

In the left panel of Fig. [15] we present an emission line Ha map, from which a properly scaled 
R band image, obtained with the same camera, had been subtracted. The morphology observed in 
the [SII] line map is very similar, showing the same t wo major features aroun d regions B and C. 
The [SII] emission at region C was already reported bv lEisloffel &: Mundtl (jl997l ) during a quiescent 
period of V1647 Ori, and it is very likely associated to the object HH22. We detected only knot A 
in our [SII] images, the other knots B to E are not seen, very likely because our images are not deep 
enough. The emission at ~ 10^ ^ North of V1647 Ori (region B) was not present in the [SII] map 
presented bv lEisloffel Mundtl (|l997l ). and it can be associated with the outburst of our target. 



We propose that the gas responsible for the emission of both Ha and [SII] is ionized by the 
hot radiation coming from the inner parts of the outburst region. In this scenario the density of 
the gas should not exceed the critical density of [SII] to avoid collisional de-excitation. Pure light 
reflection can be excluded, partly because the line emitting area is less extended than the region 
observed in the adjacent continuum (e.g. Rc or Ic band, see Fig. fTS]) . and because optical spectra 
of the central part revealed Ha but not [SII] emission. 

In the right panel of Fig. [15] we display a region located North of McNeil's nebula. The 
compact object se e n at ADec ~ 160" corresponds to knot A of the object HH23, according to 
Eisloffel fc Mundtl (Il997l ). The object is clearly seen in the [SII] map, but is faint in Ha emission, 
indicative of shock ex citation in this region. Th e source V1647 Ori was considered as the source of 
the HH23 outflow bv lEisloffel fc Mundtl (|l997l ^. 



5.3. Polarization map 

Near infrared polarization measurements were obtained using LIRIS in the J band. The po- 
larization vectors (their orientation and degree of polarization) are presented in Fig. [16] 

The general polarization pattern is centrosymmetric around V1647 Ori. The largest values 
of the polarization (slightly above 15%) are observed along the North-East rim of the nebula and 
also, remarkably, to the South- West of the star. This distribution is consistent with dust scattering 
taking place in an outer envelope surrounding and obscuring the star. Region B, the NW plume 
close to the nucleus, shows no significant polarization. There is also a conspicuous lane of low 
polarization crossing the star from East to North-West, oriented almost perpendicularly to the 
direction of region C. 



This pattern of polarization has been observed in other bipolar or cometary nebulae around 



- 29 - 



young stellar objects (jTamura et al. Ill99ll ). It is u sually attributed to the dusty dis k surrounding 
the star and it has been successfully simulated by iFischer. Henning &: Yorkd (|l996l ). Comparison 
with their simulations suggests that our observations are better explained with models which assume 
dust particles of 0.1 to 1 ^m size. Among the different accretion disk types our observations seem 
to be better reproduced by their SL model, which i nclud es a massive self-gravitating disk, in the 
form of an infinite isothermal slab (iLazareff et al.l 119901 ). The polarization towards V1647 Ori 
(3.2 % at PA = 140°) is also c onsistent both in value and orientation with the simulations by 



Fischer. Henning &: Yorkd (|l996l ) for the same accretion disk type for the unresolved polarization. 



5.4. The very red object near HH 22 



Two bright red and compact sources appear in our Kg image (Fig. [3]) towards North-East of 
V1647 Ori. The two sources are roughly aligned with the star V1647 Ori. They are detected in 
all our Ks images obtained with LIRIS, even during the post-outburst phase when the brightness 
of V1647 Ori has decreased by more than 2 magnitudes. These sources are not detected in our 
TCS images due to the poorer sensitive of the camera. The closer source is located at the position 
(Aq, A(^) ~ (8", 19"). It was not detected in our J images and only barely detected in the H 
images. Its estimated magnitudes are J > 19, > 16, and Kg ~ 15.0. Given the uncertainties in 
the available data, we cannot say more about the nature of this object. 

The farther source is located at an offset of (Aa, A(^) ~ (16", 38"). We estimate the following 
magnitudes for this source: J > 19, -ff > 16, and Kg ^ 14.5. We have checked the photometry 
of the source at different epochs, showing a variation smaller than 0.3 mags. This source has a 
corresponding entry in the 2MASS catalog (J05461423-0005261), whose measured magnitude is 
similar to ours (AKs<j — 0.4), given the uncertainties inherent to the faintness of the source. Its 
near infrare d colors (H — Kg ~ 1.5 and J — H > 3) are close to the locus of the classical T 



locus. Muzerolle et al 



Tauri stars (IMever et al.lll997l ) , although some reddening is also needed to bring its colors onto this 



(|2004l ) detected an infrared source at 4.5 and 24 /im using the Spitzer Space 
Telescope. The position of this source, as overplotted in Fig. O shows a perfect coincidence with 
the bright source in our Kg image. In the b est seeing images this source appears sligthly extended, 



as already noted bv lMuzerolle et al.l (j2004l ) 



Muzerolle et al.l (120041 ) claimed that this red object is part of the Herbig-Haro flow HH 22, and 



they propose to classify it as an extreme embedded cl ass I source. In Fig. [15] we overplotted the 
position of the red source using the coordinates given in lMuzerolle et al.l (|2004l ) on our Ha emission 
line map. It can be seen that the red source is located 3-4" North of the H« em ission line region 
which coincides with the position of knot A in HH22 (jEisloffel &: Mundtl 119971 ). Thus if the red 
Spitzer source is considered to be the driving source of HH22, it should be explained why it is not 
aligned with the chain of the knots. 



According to 



Briceho et al 



(120041 ) the location of their region C is almost centered on the 
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position of HH22. The fact that region C brightened during the outburst could be explained by 
the overlap of two regions: a dust cloud which reflects the brightening of the star, and the other 
related to HH22 producing the line emission. 



6. Discussion 



The eruptions of young stellar objects are commonly explained as a consequence of dramat- 
ically enhanced accretion from the circums tellar disk onto th e star. Based on this idea, detailed 
models for FUors have been developed bv iBell &: LinI (| 19941 ) (hereafter BL m odel), which were 
later extended to include trigger mechanisrn for the eruption (jBell et al.lll995l ). and to take into 
account reprocessed light (jTurner et al.lll997l ). In this section we first discuss whether the different 
types of observational results on V1647 Ori, collected during the two years of the outburst, could 
be consistently interpreted in the framework of the BL model. Then we comment on the question 
if V1647 Ori could be considered to be a member of either the FUor or the EXor group. 



6.1. The outburst mechanism 



Pre-outburst phase. In the BL model FUor eruptions occur in low mass pre-main sequence 
stars at an early evolutionary phase {tage < lO^yrs). These young objects are still surrounded by 
an infalling envelope which feeds material onto an accretion disk. V1647 Ori seems to be precisely 
this type of star: based on pre-outburst measurements, it is a Class I/Class II obje ct {tage ~ lO'^ yr), 
whose flat SED shows that the object is deeply embedded (jAbraham et al.ll2004l ). 



The rate at which material spirals within the disk towards the star is a critical input parameter 
of the BL model. If accretion exceeds a threshold value of about 5 x 10^'^MQ/yr, then most of 
the matter would not fall directly onto the star due to the inefficient outward transportation of 
angular momentum. It piles up in the inner disk, soon or later leading to an outburst. In Sect.H] 
we determined the accretion rate of V1647 Ori from the Br7 line fluxes for May 2006 when the star 
was in quiescent phase , and obtained a value of 5 x 10~^Mq /yr (a similar figure was proposed by 
Muzerolle et al.l (|2004l )). The result agrees with the threshold value, indicating that V1647 0ri just 
passed the criterion of the BL model to produce an outburst, but its accretion rate in the disk was 
probably the lowest among all known FUors. 



The outburst. According to the BL model, the inflowing matter piles up in the inner disk at 
^«imit~0.25 AU, until its column density - and thus opacity - becomes high enough to switch on a 
thermal instability. An ionisation front propagates in the inner disk, and the emission of the bright 
ionised region behind the front causes the increase of flux at optical and near-infrared wavelengths. 
V1647 Ori brightened by about 4.5 magnitudes in the Ic band, which is a typical value for FUor 
eruptions. One should note, however, that the brightening of V1647 Ori consisted of two effects 
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of comparable amplitude: an intrinsic brightening presumably related to the appearance of a new 
hot component in the system (probably the ionised inner disk): and a dust-clearing event which 
reduced the extinction along the line-of-sight (jReipurth &: Aspin II2004I : iMcGehee et al.ll2004l ) . Thus 
the amplitude of the outburst of V1647 Ori was probably lower than typical of FUors. Also, the 
bolometric luminosity increased only by a factor of 15, in contrast with the factor of ~100 of classical 
FUors. Concerning the hypothesis of the BL model that matter piles up at Rumit, it is interesting 
to recall that in Sect. [3] we speculated that the periodic optical dimming could be the consequence 
of variable extinction due to obscuration by a dense circumstellar dust clump orbiting the star. The 
period of the optical variability of 56 days implies an orbital radius of R = 0.28 (M/Mq)^/^ AU, 
which is close to the expected Rumit of the inner disk where material piles up, suggesting that the 
extinction might be related to this reservoir of matter. 

During the outburst the high temperature in the ionized part of inner disk leads to higher 
efficiency of angular mornentum transport, resulting in dramatically increased accretion rate onto 
the star. iMuzerolle et al.l (|2004l ) determined that the accretion rate increased by a factor of 15. In 
Sec. m we calculated an accretion rate of 5 x 10~^Mq / yr fro m th e Br7 line lumino sity for the peak 



of the outburst (close to the estimates of lVacca et al.l (j2004l ) and lGibb et al.l (j2006l )). These results 
suggest that in the eruption of V1647 Ori both the peak value of the accretion rate and its increase 
with respect to the quiescent phase were relatively modest compared to those of other FUors. 

Many observations indicated that a new, hot component appeared in the V1647 Ori system 
during the outburst. This component is probably the inner ionised part of the disk, for which the 
BL model predicts a surface temperature of 6000-8000 K. The observations seem to be consistent 
with this expectation : iMcGehee et al. I booi ) attempted to fit the SED of this new component with 
a 7000 K blackbody; IMuzerolle et al.l (|2004l ) assun ied a temperatu r e of 6 000 K for the innermost, 
outbursting part of the disk; and the disk model of [Abraham et al.l ( 20061 ) predicted a temperature 
of 4200 K for the same part. The high temperature of the inner disk is also supported by the 



CO excitation d iagrams (2400 K, iRettig et al. II2005I ). It is also consistent with the model of 



Bell fc Lid (119941 ) , that the region participating in t he outburst is liinited t o within a few AU around 
the star (see the i nterfe rometric measurements by Abraham et al. ( 2006 ): and the estimations of 
Tsukagoshi et al.l (|200,5l )). 



The main appar ent contradi c tion t o the BL model is related to the timescales. In their standard 
self-regulated model, iBell &: LinI (| 19941 ) computed a rise time of several decades (matching the case 
of the FUor V1515Cyg). In order to explain shorter initial brightening curves (e.g. V1057Cyg, 
FU Ori), they proposed 'triggered eruption' when a nearby compani on can temporarily i ncrease the 
accretion rate. V1647 0ri went into outburst in about 4 months (jBriceho et al.ll2004l ). thus it is 
tempting to speculate on such a trigger mechanism; but so far there is no observational signature 
of any companion. One can realize, however, that in the BL model all timescales depend on the 
accretion rate in the outer disk: lowering the accretion rate leads to shorter eruptions. Thus the 
very low accretion rate of V1647 Ori (practically at the threshold value for producing a n outburst) 
can be one factor to explain the short initial brightening. Table 2 in the paper of iBell &: Lin 
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(j 19941 ) presents rise time estimates for several Mj„ values. Interpolating within this table suggests, 
however, that the resulting timescales are still too long for V1647 Ori. The other main factor which 
determines the timescale of the initial brightening is the viscosity paramet er in the ou t burst high 
state of the inner disk ah- Guessing again from the numbers of Table 2 of iBell &: LinI (|1994| ). we 
concluded that a value of a/j = 3 x 10~^ (rather than their standard value of a/j = 10~^) would 
produce the observed rising timescale of about 4 months. 



The plateau phase. The initial brightening of V1647 Ori was followed by a 2-year long slow fad- 
ing. This is the period under which the inner disk depletes and the matter is becoming neutral again. 
The length of this period is determined by the same viscosity parameter as we discussed in relation 
to the length of the initial brightening. It is encouraging, that adopting again = 3 X 10-2 and 



taking into account the low accretion rate of V1647 Ori, from the numbers in Tab. 2 of lBell &: Lin 



(| 19941 ) we estimated a plateau phase of about two years, consistently with the observations. 



According to the model of iBell &: LinI (| 19941 ). the fading is due to the decreasing accretion 



rate. Checking two independent indicators of the accretion rate, in Seed] we found the somewhat 
contradictory result that the temporal evolution of the flux of the HI lines indeed shows such a 
monotonic drop in the accretion rate, while the Call 8542 A line fluxes remained constant during 
the whole outburst. 

The fading rate during this plateau was about 0.04 mag/yr. As mentioned already in Sect. 3.1, 
the fact that the fading rate appears wavelength independent suggests that no significant temper- 
ature variation occurred in the inner parts of the star-disk system in this period. 



Final fading. In October 2005, V1647 0ri started a sudden fading which took approximately 4 
months and ended when the object returned to the quiescent, pre-outburst state at optical and 
near-IR wavelengths. Remarkably, the rate of the brightness change was similar during the final 
fading and the initial brightening. This is a natural fact in the BL model, where the ionization front 
reverses at the end of the outburst and proceeds with the same speed as in the extension phase. 
We note that such a return of a FUor to the quiescent phase has never been observed before. 



Recurrent outbursts? Following an eruption, steady accretion from the outer disk slowly fills 
up the inner disk at Rumit, and the system is prepared for a new outburst. The timescale between 
outbursts is determined by the viscosity parameter of the model i n the quiescent cold state of 
the system ac- From a rough interpolation in Tab. 2 of iBell &: LinI (| 19941 ) one may conclude that 
Uc ~ 1.5 X lO"'^ would produce the 40 yr period between the last and the present outburst of 
V1647 Ori. This number is larger than the corresponding value in the standard BL model by a 
factor of 15. 
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Dust sublimation/condensation. Similarly to the initial brightening, the final fading also 
consisted of two events: an intrinsic fading and an increasing extinction along the line-of-sight. 
This suggests that the effects changing the extinction must be reversible. This reversibility is als o 
supported by the fact that V1647 0ri had already an outburst in 1966-1967 (jAspin et al.ll2006l ). 
Thus among the possible mechanisms responsible for the dust-clearing event sublimation of dust 
particles at the beginning of the outburst is a likely explanation. Is could be a reversible process, 
since the increasing extinction observed during the final fading can be due to the condensation of 
dust grains (dust condensa t ion eff ects were already observed in young eruptive stars, e.g. V1515 Cyg 
in 1980, see lKenvon et all (|l99ll )l. 



Wind. The analysis of spectral lines (Sect.|l|) suggests the presence of a wind which was strongest 
at the peak of the outburst. During the plateau phase the P Cygni profile of the prominent spectral 
lines disappeared, also indicating that the strong wind pres ent at the peak of the outburst decreased 
considerably. Such a wind is not included in the model of lBell &: LinI (| 19941 ). but should be taken 
into account for a detailed modelling of V1647 Ori. The detailed model should also take into 
account the interesting finding that the Hel absorption turned into emission in the quiescent phase, 
indicating that the type of the wind changed after the end of the outburst. 



Conc lusions. From the comparison of observations of V1647 0ri with the model of iBell &: Lin 



(| 19941 ) we conclude that the BL model reproduces most observational results, but in order to match 
the timescales one order of magnitude higher viscosity parameters ~ both in the cold and the high 
states - have to be assumed. The physical reason behind the increased viscosity parameter cannot 
be deduced from the available data. 



6.2. FUor or EXor? 



In the literature there is an on-going debate on the classification of V1647 Ori, presenting a 
number of arguments supporting either the FUor- or the EXor-like nature. On the one hand, 
in our previous paper ( Abraham et al. 2004 ) we argued that the shape of the SED is similar to 
FUors, and now we can add that the model of lBell &: LinI (119941 ). developed for explaining the FUor 
phenomenon, can model the outburst of V1647 Ori reasonably well. On the other hand, the short 
timescale of the eruption (2 years) and the recurrent nature argue more for an EXor-type event. 
EXors, however, are assumed to be classical T Tauri stars, while V1647 Ori is obviously in an earlier 
evolutionary phase than T Tau stars. We sp eculate that yi647 Ori - together with another recently 
erupted young stellar object OOSerpentis (jKospal et al.ll2006l ) - might form a new class of young 
eruptive stars. Members of this class may be defined by their relatively short timescales, recurrent 
outbursts, modest increase in bolometric luminosity and accretion rate, and an evolutionary state 
earlier than that of typical EXors. 
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7. Summary 

Comparison of our optical and near-infrared data obtained in 2004 February-2006 Sep on 
V1647 Ori with published results led to the following new results: (1) The brightness of V1647 Ori 
stayed more than 4 mag above the pre-outburst level until October 2005 when it started a rapid 
fading. (2) In the high state we found a periodic component in the optical light curves with a period 
of 56 days. (3) The time delay between the brightness variations of the star and a nebular position 
corresponds to an angle of 61° ± 14° between the axis of the nebula and the line of sight. (4) The 
overall appearance of the infrared and optical spectra did not change in the period 2004 March- 
2005 March, though steady decrease of HI emission line fluxes could be observed. We show that 
the periodic variations of the Ha equivalent width, suggested by several authors, could be caused 
by variation of the continuum with the 56 days period. In 2006 May, in the quiescent phase, the 
Hel 1.083 nm line was in emission, contrary to its deep blueshifted absorption observed during the 
outburst. (5) The J — H and H — Kg color maps of the infrared nebula reveal an envelope around 
the star whose largest extension is about 18"(0.03pc). (6) The color distribution of the infrared 
nebula suggests reddening of the scattered light inside a thick circumstellar disk. (7) Comparison 
of the Kg and Ha images of McNeil's Nebula shows that HH22A, the Spitzer infrared source 
and the bright clump C of the nebula may be unrelated objects. (8) We show that the observed 
properties of V1 647 Ori could be interpreted in the framework of the thermal instability models of 
Bell et al. (jlQQsl ). (9) We speculate that V1647 Ori might belong to a new class of young eruptive 



stars, defined by relatively short timescales, recurrent outbursts, modest increase in bolometric 
luminosity and accretion rate, and an evolutionary state earlier than that of typical EXors. 
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A. The 2MASS source J05461162 0006279 



The object located at about 30" West and 30" S of V1647 Ori (Fig.ED can be identified with 
the 2MASS source J05461162— 0006279. The source is very re d as shown by its app earance as a 
remarkably bright object in Spitzer im ages at 4.5 and 24 fim (jMuzerolle et al.ll2004l ). It was also 
detected by XMM (iGrosso et al.ll2005l ). however it remains undetected in optical images. Thanks 
to the excellent seeing (0'.'65) during the night 2004 Nov 4 we could resolve this target into two 
components separated by 1"2, oriented along the E-W axis (see Fig. The magnitude of each 
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component was obtained by performing PSF photometry using the daophot package in IRAF. The 
primary star is the Western component, and has J= 14.14 and = 11.43; and the secondary star 
has J= 15.72 and = 13.79. The color index J— is 2.7 and 1.9 for the primary and secondary 
stars, respectiv e ly. T hese color indices are redder than those of normal M stars according to 



Bessell fc: BrettI (j 19881 ). suggesting the possibility that this double system is also embedded in the 



Orion B molecular cloud. 
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Fig. 11. — Near-infrared {HK band) spectra of V1647 Ori obtained with LIRIS, normalized to the 
continuum. 
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Fig. 12. — Flux evolution of near-infrared and optical spectral lines of V1647 0ri. Circles: this 
study, open circle s correspond to Hel emission component (Tab.Bl) ; dia monds: dat a take n from 



Gibb et al.l (|2006l ): triangles: data taken from IWalter et al.l (|2004l ) and lOjlia et al.l (|2006l ). The 



formal error bars are smaller than the symbol sizes. 
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Fig. 13. — Color maps showing the V — Rc {left panel) and H — Kg {right panel) indices of the 
infrared nebula around V1647 Ori. The contours in the left and right panels represent the emission 
in Rc and J bands, respectively. 




Fig. 14. — Different color profiles of the nebula along N— S, crossing the star. In the left panel 
we represent the variation of the V — Rc (stars) and R — Ic (crosses) colors, underlying is the 
brightness profile as measured in the R band. In the right panel we represent the variation of the 
J — H (diamonds) and H — Kg (triangles) colors, underlying is the brightness profile as measured 
in the J band. 
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Fig. 15. — Left panel: The color scale image corresponds to the Hq emission after continuum 
subtraction. The contour lines represent the emission in the Rc continuum. Th e filled diamond- 



marks the position of the bright infrared source dete cted in our K,^ b and a nd also bv lMuzerolle et al 



(j2004l ) in Spitzer images. Apertures B and C from lBriceho et al.l (|2004l ) are also indicated. Right 
panel: The color scale image corresponds to the Ha emission in the region around the Herbig-Haro 
object HH 23. The contour lines represent the [SII] emission. Note the bright knot detected in [SII] 
emission at ~ 160" N, which appears very faint in the Ha. The scale is twice that of the left panel. 
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Fig. 16. — Polarization map obtained in the J band. The color scaled image represents the total 
surface brightness distribution. 



